Maple 6 Questions and Posts

These are Posts and Questions associated with the product, Maple 6

Can anyone help me?

restart: with(plots):   

h:=0.01: n:=6000: Digits:=18:

t[0]:=0: x[0]:=0.2: y[0]:=0.1: epsilon:=20:

F:=(x,y)->epsilon*(1-x^2)*y-x:

begin:=time():

for k from 0 to n do

t[k+1]:=t[k]+h;

x[k+1]:=x[k]+h*y[k];

y[k+1]:=y[k]+h*F(x[k],y[k]);

pt[k]:=[t[k],x[k],y[k]]:

end do:

cpu_time:=time()-begin;

pointplot3d([seq(pt[j],j=0..n)],axes=normal,symbol=cross,symbolsize=8,color=red,labels=["t","x","y"],orientation=[-90,0]);

sys:=diff(X(t),t)=Y(t),diff(Y(t),t)=F(X(t),Y(t));

vars:={X(t),Y(t)}: ic:=X(0)=0.2,Y(0)=0.1:

sol:=dsolve({sys,ic},vars,numeric,stepsize=h,method=classical[foreuler],output=listprocedure):

odeplot(sol,[t,X(t)],0..h*n,axes=normal,style=line,numpoints=n,labels=["t","x"]);

XX:=eval(X(t),sol): x1:=XX(1);  pt[100];

 

 

Download 09-1-1.mws

Hello!

I have a function, say f(x,y), and I want to solve for x for a range of values of y. The function is complicated, and thus, I gather I need to use fsolve to find the set of solutions numerically. I know how to use fsolve to find x for a given (single) value of y, but I have never used it to generate a list of x values given a range of y values. Then, I would l like to plot x as function of y, but I will leave that for later.

 

Thank you!

I have been using fsolve to derive constants for a polynomial. These are presented clearly in blue, but I am unable to reference these numbers for further calculations other than by typing them in manually. e.g. output is:

{d = 520.9032305, c = .41478630e-1, b = -.84548378e-5, a = .125070852e-8}

 


fsolve({y[1]=a*x[1]^3+b*x[1]^2+c*x[1]+d,y[2]=a*x[2]^3+b*x[2]^2+c*x[2]+d,
> y[3]=a*x[3]^3+b*x[3]^2+c*x[3]+d,y[4]=a*x[4]^3+b*x[4]^2+c*x[4]+d},{a,b,c,d});

can anyone one help me to download maple 6 ?

I've used fsolve in Maple 6 to find the roots (real and complex) of a polynomial.  I'd like to create a plot showing the location of those roots in the complex plane.  How can I do this?

Page 1 of 1